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ABSTRACT 

Radio properties of supernova outbursts remain poorly understood despite 
longstanding campaigns following events discovered at other wavelengths. After 
~ 30 years of observations, only ~ 50 supernovae have been detected at radio 
wavelengths, none of which are Type la. Even the most radio-loud events are 
~ 10 4 fainter in the radio than in the optical; to date, such intrinsically dim 
objects have only been visible in the very local universe. The detection and 
JpH| study of radio supernovae (RSNe) will be fundamentally altered and dramati- 

cally improved as the next generation of radio telescopes comes online, including 
EVLA, ASKAP, and MeerKAT, and culminating in the Square Kilometer Array 



(N 



(SKA); the latter should be > 50 times more sensitive than present facilities. 
SKA can repeatedly scan large (> 1 deg 2 ) areas of the sky, and thus will discover 



>- ' RSNe and other transient sources in a new, automatic, untargeted, and unbiased 

way. We estimate SKA will be able to detect core-collapse RSNe out to red- 
shift z ~ 5, with an all-redshift rate ~ 620 events yr _1 deg -2 , assuming a survey 
sensitivity of 50 nJy and radio lightcurves like those of SN 1993J. Hence SKA 
should provide a complete core-collapse RSN sample that is sufficient for statis- 
tical studies of radio properties of core-collapse supernovae. EVLA should find 
~ 160 events yr _1 deg -2 out to redshift z ~ 3, and other SKA precursors should 
have similar detection rates. We also provided recommendations of the survey 
strategy to maximize the RSN detections of SKA. This new radio core-collapse 
supernovae sample will complement the detections from the optical searches, such 
as the LSST, and together provide crucial information on massive star evolution, 
supernova physics, and the circumstellar medium, out to high redshift. Addi- 
tionally, SKA may yield the first radio Type la detection via follow-up of nearby 
events discovered at other wavelengths. 
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Introduction 



Supernovae are among the most energetic phenomena in the universe, and are central 
to cosmology and astrophysics. For example, core-collapse supernovae are explosions arise 
from the death of massive stars and hence are closely related to the cosmic star-formation 
rate and to massive- star evolution, and are responsible for the energy and baryonic feedback 
of the environment (IMadau et al.lll998l ). Type la supernovae show a unifor m properties in 
their lightcurves and play a crucial role as cosmic "standardizable candles" ( jPhillipsI Il993l ; 
Riess et al.lE993 ). 



Our knowledge of the optical properties of supernovae, is increasing rapid 
vent o f prototype "synq ptic"-i.e., rep eated scan-sky surveys, such as SDSS-II ( 



20081 : ISako et all 12008 1 and SNLS (IBazin et al.l 120091 ; IPalanque-Delabrouille et al.1 BOlOh . 



y with the ad- 



Frieman et al. 



These campaigns are precursors to the coming "Great Survey" era i n which synoptic surveys 
will be conducted routinely over very large regio ns of sky, e.g., LSST ( The LSST Collaboration 
20071 ILSST Science Collaborations et alil2009l ) and Pan-STARRS fabriry] 120031) The num- 



ber o f detected supernovae will increase by several orders of magnitude in this era (ILSST Science Collaborate 
2009J ; JLien fc Fields! l2009h . 



In contrast to this wealth of optical information, properties of supernovae in the radio 
remain poorly understood, fundamentally due to observational limitations. Radio super- 
novae (RSNe) have primarily been discovered by follow-up observations of optical outbursts, 
and only very rarely by accident. To date, only ~ 50 core-collapse ou tbursts have radio 



detections, and no T ype la explosion has ever been detected in the radio (IWeiler et al.ll2004 



Panagia et al.ll2006l ). The core-collapse subtype Ibc has been a focus of recent study in the 
radio, because some Type Ibc events are associated with long Gamma- Ray B ursts (GRBs) 
jGalama et alill998[ JKulkarni et aJ]ll998[ ISoderbereJl2007l ; berger et aliboOoJ 



Current radio interferometers are scheduled primarily around targeted observations pro- 
posed by individual principal investigators. This stands in contrast to future radio interfer- 
ometers planned for the coming "Great Survey" era. These include the Square Kilometer 
Array (SKAlj) and its precursor prototype arrays (for example, ASKAPg and MeerK ATjj) . 
These telescopes will operate primarily as wid e- field survey instruments focusing on several 
key science projects (jCarilli fc Rawlingsj 120041 ) . As synoptic telescopes, they will be far bet- 
ter suited to study all classes of transient and time-variable radio sources, including RSNe. 



1 http://www. ska-telescope. org 
2 http://www. atnf.csiro.au/projc 
3 http://www.ska.ac.za/meerkat 
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Gal-Yam et al.l ( 120061 ) already pointed out the power of synoptic radio surveys for detecting 
radio transients of various types, including supernovae and GRBs, in an unbiased way. Here 
we quantify the prospects for RSNe. 

In this paper we explore this fundamentally new mode of untargeted RSN discovery and 
study. We adopt a forward-looking perspective, and consider the new science enabled by 
RSNe observations in an era in which the full SKA is operational. Our focus is mainly on 
core-collapse supernovae, the type for which some radio detections already exist. However, 
we will also discuss the possibility of Type la radio discovery based on current detection 
limits. We will first summarize current knowledge of radio core-collapse supernovae (§2}, 
and the expected sensitivity of SKA (§S\). Using this information, we forecast the radio 
core-collapse supernovae harvest of SKA (§1J, and consider optimal survey strategies (§5]). 
We conclude by anticipating the RSN science payoff in this new era (§[§]). We adopt a 
standard flat ACDM model with fi m = 0.274, fi A = 0.726, and H = 70.5 km s" 1 Mpc -1 



( JKomatsu et al.l 120091 ) throughout. 



2. Radio Properties of Supernovae 



Several key properties of RSNe have been establish ed, as a result of the longstanding 



leadership of the NRL-ST ScI group (recently reviewed in lWeiler et al.ll2009 



Stockdale et al. 



2007 : Panagia et al.ll2006l) and of the CfA group and others (summarized in lSoderbergl 



2007 



Berger et al.l 120031 ) . We summarize these general RSN characteristics, which we will use to 
forecast the RSN discovery potential of synoptic radio surveys. 



2.1. Radio Core-Collapse Supernovae 



Observed core-collapse RSNe have luminosities spanning vL v ~ 10 33 — 10 38 erg •sT 1 at 5 
GHz, and thus are > 10 4 times less luminous in the radio than in the optical. Their intrinsic 
faintness has prevented RSN detection in all but the most local universe. Even within a 
particular core-collapse subtype, radio luminosities and lightcurves are highly diverse, e.g., 
two optica lly similar Type Ic events might be radio bright in one case and undetectable in 
the other dMunari et al.lll998[ iNakano fc Aok3ll997l : Istockdale et~aD J2006h 0. Additionally, 
core-collapse RSNe spectral shapes strongly evolve with time; lightcurves peak over days to 



4 New Radio Supernova Results (|Stockdale et al. 



20061) 



are available online at: 



http://rsdwww.nrl.navy.mil/7213/weiler/sne-home.litml 
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months depending on the frequency. RSN emission can be understood in terms of in teractions 
betwee n the b last, ambient relativistic electrons, and the circumstellar medium ( [Chevalier 
1982al lbl. [l998h . 



To model R SN emission as a f unction of frequen cy and time, we ado pt the semi-empirical 
form derived by I Chevalier! (jl982al ) and extended in IWeiler et al.l (120021 ) , 



L(t, v) = L x 



5 GHz 



1 day 



'external 



e -TCSM clumps 



r CSM 



clumps 



^internal 



'internal 



We follow the notation of IWeiler et al.l (J2002h . L(t,u) is the supernova luminosity at fre- 
quency v and time t after the explosion. Optical depths from material both out side (r PYternfl .i , 
tcsm t ) and inside (Ti ntcrna i) the blast-wave front are taken into account (see IWeiler et al. 
20021 ). 



Parameters em bedded in each op tical depth term are those for SN 1993J, one of the 
best studied RSNe (IWeiler et al.ll2007l ). Radio emission from SN 1993J is dominated by the 
clumped-circumstellar-medium (clump-CSM) term, and hence 

1 p~ TCSM clumps 

Lfau) , (2) 

'CSM clumps 

T ^)- 2 - 83 , for SN 1993J. At small t, r CS M clumps is 



where tcsm , 
large and L(t, v) 



4.6 x 10 E 



,-2.1 



-5 GHz' 



l/ r csM clum s oc v 2 - 1 t 2 - 83 , so luminosity grows as a power law at early 



times. With all optical depth parameters fit to SN 1993J, the peak luminosity is controlled 
by the prefactor L\. 

Our main focus will be on RSN discovery, and thus it is most important to capture the 
wide variety of peak radio luminosities, which correspond in our model to a broad distri- 
bution for L\. Figure [1] shows a crude luminosity function (not-normalized) based on the 



sample of 20 core-collapse super novae ( 



peak luminosity at 5 GHz (IWeiler et al. 



5 Type II and 5 Type Ibc ) that have a published 



2004 : IStockdale et all 1200.1 20071 : IPapenkova et al. 



200ll : IStockdaleet al.l 120061 : iBaklanov et al.ll2005l : IPoolev et al.ll2002h . We use 5 GHz data 
to construct our luminosity function because the most RSNe have been observed at this 
frequency. However, our predictions will span a range of frequencies, based on this lu- 
minosity function and eq. (J2j) . The data are divided into four luminosity bins of size 
Alog 10 (L) = 1. The black curve in Fig. [T] is the best-fit Gaussian, with average lumi- 
nosity log 10 (L avg /ergs _1 Hz -1 ) = 27.3, a standard deviation a = 1.25, and \ 2 — 0.18. SN 
1987A is marked in Fig. [H but was not used in the fit to avoid possible bias due to its 
uncommonly low luminosity. The fitted luminosity function might be biased towards the 
brighter end, because of the current survey sensitivity and the small and incomplete nature 
of the sample. 
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Fig. 1. — Radio luminosity function (not-normalized) at 5 GHz of core-collapse supernovae 
showing core-collapse supernovae count as a function of log 10 (L), where L is the peak radio 
luminosity. Data are binned to A log 10 (L) = 1. The black curve shows the x 2 -htted Gaussian 
to the underlying data (red stars). 



2.2. Radio Type la Supernovae 

All searches to date have failed to detect radio emission from Type la supernovae. 
Panagia et al.l (120061 ) reported the radio upper limits of 27 Type la supernovae from more 
than two decades of observations by the Very Large Array (VLA). The weakest limit on a 
Type la event is 4.2 x 10 26 erg s _1 Hz -1 at 1.5 GHz for SN 1987N, which is around one order of 
magnitude lower than the average luminosity of radio-detected core-collapse supernovae (see 
§2.ip . The strongest limit on Type la radio emission is even tighter, 8.1 x 10 24 erg s _1 Hz -1 
at 8.3 GHz for SN 1989B. Additionally, t he z ~ cosmic T ype la supernova rate is around 
1/4.5 of the core-collapse supernova rate ( IBazin et al.ll2009[ ). The intrinsic faintness in radio 
and their smaller rate make detecting Type la in radio observations especially hard. 



3. Next- Generation Radio Telescopes: Expected Sensitivity 

Radio detections of supernovae to date have been restricted both by the limiting sensi- 
tivity of contemporary radio interferometers and the need for dedicated telescope time for 
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transient followup. This situation will change drastically with SKA's unprecedented sensi- 
tivity and particularly by its ability to repeatedly scan large regions of the sky at this great 
depth. 

Current SKA specifications adopt a target sensitivity parameter A e s/T sys = 10 4 m 2 K _1 
at observing frequencies in the low several GHz, including z = HI at 1.4 GHz. A e s is the 
effective aperture, and T sys is the system temperature. We will adopt this value of v4 e ff /T sys , 
which yields a 1-a rms thermal noise limit in total intensity of 

oi = 0.15 /iJy (Azz/GHz)- 1 / 2 {5t/hr)- 1/2 , (3) 

for a bandwidth Av and observation duration St. The SKA will therefore reach a thermal 



noise limit of several nJy in deep continuum integrations (St ~ 1000 hr) (ISKA Design Reference Mission 



20091 ) [3. We define the associated survey sensitivity S m - m (the minimum flux density thresh- 



old) as ^min = 3a 1. In common with other radio interferometers, SKA will accumulate 
sensitivity in targeted deep fields, including transient-monitoring fields, by accumulating in- 
tegration time over multiple individual observing tracks. We therefore will adopt a fiducial 
SKA supernova sensitivity of S m { n = 50 nJy in 100 hours of observation, but we will show 
how our results are sensitive to other choices of S m i n . 

It is anticipated that transient fields will be revisited with a cadence appropriate to 
the variability timescales under study and that interferometric inverse imaging methods will 
include source models with time variability. Survey optimization for interferometric transient 
detection is an active area of current SKA research. The technical details are beyond the 
intent and scope of this paper, but will be influenced by science goals for transient source 
study in general, including pulsars, GRBs, and supernovae (as considered in this paper), as 
well as the as-yet undiscovered transient population. 



4. Radio Supernovae for SKA 

With its unprecedented sensitivity, SKA will be capable of synoptic search for core- 
collapse RSNe and open new possibilities in radio astronomy. In this section, we predict the 
RSN detections of SKA based on current knowledge to demonstrate how the RSN survey 
can be done. 



5 http://www.skatelescope.org/PDF/DRM_vl.0.pdf 
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4.1. Core-Collapse Supernovae 

The detection rate Tdetect = dNgN / (dt ^ s dfi dz) for a given RSN survey is 



J- detect /survey /radio /ISM J- 



SN 



(4) 



and is set by several observability factors / that modulate the total rate of all supernovae 



Tsn(^) 



dN< 



SN 



^pm Ct V rt 



& 'comov ^em ^^obs ^ ' UZ 



KsN(z)r, 



2 ( 

comov ' 



z) C 



dt 



dz 



(5) 



within the cosmic volume out to redshift z (IMadau et al.l 11998c iLien fc Fieldsl 120091 ) . 



We see that the total cosmic supernova rate Tsn depends on cosmology via the vol- 
ume element and the time dilation terms. Because ACDM cosmological parameters are 
now known to high precision, these factors have a negligible error compared to the other 
ingredients in the calculation. The other factor in Tsn is the cosmic core-collapse supernova 

Some direct measurements of this rate now 

1 



rate density TZ SN (z) = dN SN / (dV comav dt emit 

exist out to z ~ 1, but the uncertainties remain large 



2004 



Cappellaro et al.l 12005c Irlopkins k, Beacoml 



2006 



ICappellaro et al. 



Botticella et al 



1999; 



2008; 



Dahlen et al. 



Dahlen et al. 



2008 



Kistler et al 



2011b 



2008 



Horiuchi et al. 



Bazin et al. 



2009; 



Smartt et al 



2009 



Dahlen et al. 



2010 



Li et al. 



201 ll ). However, core-collapse events are short-lived, and so the cos- 



mic core-collapse rate is proportional to the cosmic star-formation rate /?*, which is much 
bett er-determined and ex tends to much higher redshifts. We thus derive 7^-sn from the re- 
cent iHoriuchi et al.l ( 120091 ) fit to the cosmic star-formation rate. The proportionalit y follows 



from the choice of initial mass function; we apply the Salpeter initial mass function (jSalpeter 
19551 ) and assume the mass range of core-collapse SNe progenitors to be 8M — 5OM ; this 
gives ft SN = (0.007 Mq 1 )^ 

Several effects reduce the total rate Tsn to the observed rate r detect in eq. fll]). Due to 
finite survey sensitivity, only a fraction f survey of events are bright enough to detect, and only 
some fraction / ra dio of supernovae will emit in the radio. We neglect interstellar extinction 
and assume /ism ~ 1 at the radio wavelengths considered. 

The term /radio in eq. (jl|) contains the greatest uncertainty due to the relatively small 
sample of RSNe observed to date, and the unavoidable incompleteness of the sample (K. 
Weiler, private communication 2010). The only published frac tion available is for Type Ibc 
supernovae. Using VLA for radio follow-up. iBerger et al.l ( 120031 ) suggests that / ra dio,ibc ~ 12% 
after surveying 33 optically-detected Type Ibc supernovae. For the purpose of demonstration, 
we will adopt / ra dio = 10% for the calculations presented in this paper, which we believe is 
rather conservative. 



An order-of-magnitude calculation provides a useful estimate of the expected core- 
collapse RSN rate. As discussed in ^J3j we adopt a fiducial SKA sensitivity of S min = 
50 nJy. Hence SKA will be able to detect supernovae with average radio luminosity (L ~ 
10 27 erg s _1 Hz -1 ) to a distance D L = y/L/AirSmm ~ 4 Gpc, which for a ACDM cosmology 
corresponds to z ~ 1. This will give a detectable volume of Vdetect ~ (4/3)7rD| ~ 2.85 x 
10 11 Mpc 3 . Observations show that the core-collapse supernova rate -Rsn ~ 10~ 3 yr _1 Mpc~ 3 
at z ~ 1 ( iDahlen et al.ll2008[ 120101 ). Assuming the fraction of the total core-c ollapse super- 
novae that display the adopted average radio luminosity to be / ra dio ~ 10% (IBerger et al. 
20031 ). the all-sky detection rate dN SN /dt ~ i? S N x /radio x Vdetect ~ 2.85 x 10 7 yr _1 . This 
corresponds to a areal detection rate gL/Vsn / (dt dQ) ~ 700 yr _1 deg -2 . As we now see, a 
more careful calculation confirms this estimate. 
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Fig. 2.— Estimated radio core-collapse supernova detection rate as a function of redshift 
at 1.4 GHz, assuming /radio = 10%. Predictions are shown for different survey sensitivities: 
<Smm = {10 /iJy (blue), 1 /iJy (blue), 100 nJy (blue), 50 nJy (thick red), 10 nJy (blue), 1 
nJy (blue)} from bottom to top solid curves, respectively. We adopt 50 nJy as our bench- 
mark sensitivity hereafter. For comparison, the red-dashed cur ve shows the LSST optical 
supernova detection rate per year per deg 2 ( JLien fc Fields! 120091 ). Also, the top solid curve 
(black) plots the ideal core-collapse RSN rate for comparison. 
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Fig. 3.— Estimated radio core-collapse supernova detection rate as a function of redshift 
for different frequency bands, for / ra dio = 10%, and an adopted survey sensitivity S min = 50 
nJy. 
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Fig. 4. — Core-collapse detection sensitivity to supernova radio luminosity, at 1.4 GHz, and 
for survey sensitivity S m i n = 50 nJy. (a.) Left Panel: Supernova distribution over redshift, 
for different cutoffs for the luminosity function, (b) Right Panel: Supernova distribution in 
luminosity bins, integrated over all redshifts. 
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A careful prediction involves detailed calculation of /survey (2) • The fraction j ' SU rvey (z) of 
observable radio-emitting events depends on adopted survey sensitivity, and on the normal- 
ized supernova luminosity function <3> 5G Hz(logL), which is measured at a peak luminosity at 
5 GHz (derived in §[2]). In this paper we will only consider whether a supernova is detectable 
at its peak luminosity at each corresponding frequency. The peak radio luminosity should 
be re ached earlier at hig her frequencies because of preferential absorption at lower frequen- 



cies ( fWeiler et al.l 120021 ). At different redshift, the peak flux density S^ m in the observed 
frequency v can be tied to the corresponding luminosity threshold L ^ by 

-^min ( Z 'i ^emit) = /i, \ ^min (^obs), (6) 

where the luminosity distance is D L (z) = (1+z) c/H J Q 2 dz' [O m (l+ z') 3 +rt\]~ 1 ^ 2 . However, 
because the luminosity function we used is based on the peak luminosity at 5 GHz, we must 
find the corresponding luminosity threshold at this frequency by applying corrections based 
on the radio spectrum, 

rpeak _ j- peak J5GHz band f (^cm) ^cm ,„. 

ml n,5GHz mi n ^^ ^ 5peak[(1 + z>ohs] ^^ ■ I ) 

The detectable fraction resulting from survey sensitivity can therefore be estimated as 



/survey = / *5GHz (log L) dlogL. 



; ^min,5GH 



Figure. [2] plots the results of our predicted core-collapse RSN detection rate for different 
target survey sensitivities, S m i n . We adopt a benchmark frequency of 1.4 GHz because this 
will be one of the first major bands SKA deploys to observe neutral hydrogen. The related 
instantaneous field-of-view at 1.4 GHz of current SKA designs based on dish reflectors is 
approximately 1 deg 2 , which we adopt. Fig. |2J plots the ideal core-collapse supernova rate for 
comparison (assuming infinite sensitivity but / ra dio = 10%). One can see that the detection 
rate at 1 nJy is very close to the ideal rate in the universe. 

Results for our fiducial SKA flux limit S m i n = 50 nJy are highlighted in Fig. [2J At this 
sensitivity, we see that we can expect that radio supernovae will be discovered (event rates 
> 5 RSNyr -1 deg -2 ) over the enormous redshift range z ~ 0.5 to 5. The total rate of RSNe 
expected in this entire redshift range is 

— 4£(> 50 nJy) w 620 RSNe yr" 1 deg" 2 , (9) 

at ail 

in good agreement with our above order-of-magnitude estimate. This sample size is large 
enough to be statistically useful and to allow for examination of the redshift history of RSNe. 



- 11 - 



Moreover, out to z ~ 1, SKA will detect almost all cosmic RSNe in the field of view, while 
at higher redshift the detections still comprise > 10% of the underlying ideal cosmic rate. 
For comparison, we also see that LSST will detect optical supernovae out to z < 1. Thus 
SKA will be complementary to LSST as a unique tool for cosmic supernova discovery. 

Figure. [3] shows how core-collapse RSN detections vary for different observing frequen- 
cies, fixing a common survey sensitivity S min = 50 nJy and bandwidth Av = 1 GHz. Results 
show similar numbers of detections at different bands, which is caused by a relatively flat 
spectrum shape at peak luminosities. Because SKA will be able to detect core-collapse RSNe 
out to high r edshif t z ~ 5, the frequency- redshift and time-dilation effects are significant. 
Weiler et al.l ( 120021 ) noted that RSNe peak when the optical depth t ~ 1. Since the optical 



depth depends both on frequency and time with similar power index (jWeiler et al.ll2002l ). the 
frequency-redshift and time-dilation effects approximately cancel, so that a fixed observed 
frequency, the peak time is nearly redshift-independent. 

As mentioned above, our luminosity function is likely biased toward the available bright 
events in a small and incomplete sample. To explore how this bias could affect our results, 
Fig. H] shows how the detection rate with S min = 50 nJy at 1.4 GHz depends on core- 
collapse RSN luminosity. Fig. HJ^a) shows that RSN with peak luminosities greater than 
10 27 erg s _1 Hz -1 contribute all of the detections beyond redshift z ~ 1, and RSN need to 
peak brighter than 10 28 erg s _1 Hz -1 to be seen beyond z ~ 3. Fig. H^b) similarly shows that 
the all-redshift detection rate becomes substantial for explosions peaking > 10 26 erg s _1 Hz -1 . 

Type Ibc supe rnovae are of particular interest given their association with long gamma- 
ray bursts (GRBs; iGalama et al.lll998t IWoosleylll993l : iHeger et al.l 120031 ) . Fig. [5] shows our 
predictions for Type Ibc detections of SKA per year per deg 2 at 1.4 GHz with a survey sensi- 
tivity of 50 nJy. The red curve shows t he radio Type Ib c detections, assuming that Type Ibc 
represents 25% of core- collapse events (ILi et al.ll2011al ). and / ra dio,ibc = 12% with luminosity 



10 27 erg s l Hz l |j ( IBerger et al.l 120031 ) . The blue curve shows the possible detections of 



the sub-class of Type Ibc supernovae that display extreme radio emission and hence might 
be powered by central engines and related to GRBs. We assume that 0.5% of all Type Ibc 



supernovae are pow ered by central engines and have luminosities of ~ 10 erg s Hz 



( Berger et al. 



( 1 Weiler et al. 



200311 . We adopted the spectrum of SN 1998bw, which is a Type Ic supernova 
20021 ). Under these assumptions the SKA will be able to make unbiased, un- 

targeted detections of ~ 130 radio Type Ibc supernovae per year per deg 2 , and ~ 20 Type 

Ibc supernovae that might be connected to GRBs. 



6 Here we simply assume a Gaussian distribution for the luminosity function centered at the specified 
luminosity with a = \. 
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Fig. 5. — Predicted detection rate of Type Ibc supernovae as a function of redshift. In this 
plot we assume the sensitivity for SKA is S m i n = 50 nJy. The re d curve shows all of the 
radio Ibc detections, assuming / ra dio,ibc = 12% (IBerger et al.1 120031 ) . The blue curve shows 
only the detection rate for Radio Ibc with central engines, assuming 0.5% of all of the Type 
Ibc RSNe are powered by central engines. 

Finally, we turn to SKA precursors. The EVLA[j, a current leading-edge radio inter- 
ferometer operating at centimeter wavelengths, is anticipated to reach a 1-cr rms noise of 
(j i ~ 1 /iJy or less in 10 hours, while SKA is expected to reach 07 ~ 50 nJy in 10 hours. 
With data accumulated over repeated scans spanning over 1000 hours, an rms 07 ~ 5 nJy 
may be reached. In synoptic surveys, we would expect EVLA to detect core-collapse events 
at a rate ~ 160 RSNe yr _1 deg -2 over a redshift range z = 0.5 to 3 (Fig. [2]). A sample of this 
size over this redshift range will already mark a major advance in the study of cosmic RSNe, 
and further motivate the full S KA. ASKAP and Mee r KAT are expected t o have sensitivities 
comparable to that of EVLA (IJohnston et al.1 120091 ; Ide Blok et al.1 12010| ) , hence we would 
expect these to detect RSNe with similar rates and redshift reach. 



7 http://www. aoc.nrao.edu/evla 
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4.2. Type la Supernovae 

If all Type la RSNe are dimmer than the weakest limit presented in §2.2[ the expected 
SKA detection rate is essentially zero. For example, if a typical Type la has a radio luminosity 
equal to the lowest published limit, L — 8.1 x 10 24 erg s _1 Hz -1 , this can be seen with a 
sensitivity S m i n = 50 nJy out to a luminosity distance ~ 300 Mpc (z ~ 0.08). While ~ 3900 
cosmic la events should occur per year out to this distance over the entire sky, <C 1 events 
are expected in the SKA field of view. More optimistically, imagine a typical Type la radio 
luminosity is L = 10 26 erg s -1 Hz -1 , which is below L = 4.2 x 10 26 erg s _1 Hz -1 , the highest 



published limit (IPanagia et al.l 120061 ); here the luminosity distance increases to ~ 1400 Mpc 
[z ~ 0.28). In this case, we find an SKA Type la detection rate ~ 0. 5 yr" 1 deg~ 2 , base d 



on the local cosmic Type la rate derived from SDSS-II optical data (jDilday et al.ll2010l ). 
S'min = 50 nJy, and /radio = 10%Jj We see that even optimistically, we expect fewer than 
one event per SKA field-of-view per year. Even with /radio = 100%, the detection rate is still 
only ~ 5 yr _1 deg~ 2 . Therefore we conclude that SKA will make few, if any, blind detections 
of Type la supernovae. 

Targeted radio observations to follow up from nearby optical detections will probably 
be the best way to search for such events. For example, we expect 10 Type la events/year 
in the LSST sky within ~ 60 Mpc (z ~ 0.015). Type la (or core collapse!) events within 
this distance observed with S min = 50 nJy, would be detectable at luminosities L > 3.0 x 
10 23 erg s _1 Hz -1 . Amusingly, this is close to the radio luminosity of SN 1987A. 



5. Radio Survey Recommendations 

A key requirement for detecting weak radio emission from CSM-supernovae interactions 
is improved r adio interferometer sensitivity. High angular resolution - below an arcsecond 



at 1.4 GHz, (jWeiler et al.l 120041 ) - is also required to avoid natural confusion and to help 
identify supernovae against background galaxies. This is similar to the maximum EVLA 
angular resolution at 1.4 GHz. For comparison, the maximum anticipated SKA baseline 
length of 3000 km, producing angular resolution of 0.014 arcsecond at 1.4 GHz, is sufficient 
to distinguish different galaxies and also to resolve galaxies as extended sources within the 



obser vable universe with rms confusion limit of < 3 nJy at 1.4 GHz (jCarilli fc Rawlings 



200J). 

A key science goal of the SKA is to detect transient radio sources, both known (e.g. 



This also is implied by Fig. [4] which is for core-collapse events that have a higher cosmic rate. 
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pulsars, GRBs), and as-yet unknown. This requires sophisticated transient detection and 
classification algorithms very likely running commensally with other large surveys planned 
by the SKA, such as the HI spectroscopic survey and deep continuum fields. We assume 
here that SKA transient detection algorithms will encompass automated detection of RSNe . 



For example, cu rrent parameterized models (jWeiler et all Il986l . Il990l ; iMontes et al.l 11997 



Chevalierlll982al Jbl) based on available data predict patterns of spectral index evolution char- 
acteristic of supernovae in general, and supernova sub-types in particular. This information 
could be exploited for RSN detection, even potentially against a back ground of unrelate d 



source variability. Broad frequency coverage is important in this regard (JWeiler et al.ll2004j ). 



The SKA intrinsically is a high dynamic-range instrument, given the sensitivity implied 
by the large collecting area. The most demanding SKA science applications will require a 
dynamic range of 10 7 :1. The detection of faint RSNe will require a dynamic range that falls 
within that envelope. 

Although the lightcurves of RSNe show great diversity, the luminosities of core-collapse 
supernovae usually change much slower in radio than in optical. RSN lightcurves typi- 
cally evolve on timescales of weeks to years; a useful lightcurve compilation appears in 



Stockdale et al.l ( 120071 ) . Thus the minimum survey cadence (revisit periodicity) need not be 
any more frequent than this. Also, we have shown that core-collapse RSNe can be found out 
to high redshift with surveys pushing down to S min = 50 nJy. For SKA this corresponds to 
about 100 hours of exposure, in line with planned deep field exposures which are part of the 
key science. Thus, SKA as currently envisioned is well-suited to core-collapse discovery. 

On the other hand, SKA probably will not have sufficient survey sensitivity for a volu- 
metric search for Type la events, based on our current knowledge of the cosmic Type la rate 
and the upper limits in their luminosities set by the non-detection of these events. Follow-up 
observations from other wavelengths will likely be the best way to search for Type la RSNe. 

The small volume of the local universe will limit nearby untargeted SKA detections 
of low-redshift core-collapse RSNe. We estimate only ~ 2 core-collapse RSN detections 
per year per square degree within redshift z ~ 0.5 (assuming a 50 nJy sensitivity at 1.4 
GHz and / ra dio = 10%). Unless SKA has large sky coverage comparable to those of optical 
surveys, it will be hard to get statistical information from such a small sample. Therefore, 
targeted radio followup of optically-confirmed nearby supernovae will be crucial to build a 
core-collapse RSN "training set" database needed for refining automatic identification and 
classification techniques. 

With detection methods optimized based on low- z radio data for optically-identified 
events, radio surveys can then be used to independently detect core-collapse RSNe at high 
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redshift based only on their radio emission. As shown in Fig. [21 supernova searches at 
high redshift (z > 1) will largely rely on radio synoptic surveys, the inverse of the strategy 
proposed above for low-redshift domain. Surveys for core-collapse RSNe will likely not be 
synoptic all-sky surveys due to operational limitations, but will likely proceed in a limited 
set of sub-fields, visited over an hierarchical set of cadences to cover a range of time-scales for 
general transient phenomena and multiple commensal science objectives. It is also important 
to match core-collapse RSNe survey sky coverage and cadence to that used in complementary 
optical surveys. LSST will repeatedly scan the whole sky every ~ 3 days. Thus a cadence 
~ 1 week for RSNe sub-fields will be preferred for an SKA core-collapse supernova survey. 



6. Discussion and Conclusions 

SKA's capability for unbiased synoptic searches over large fields of view will revolu- 
tionize the discovery of radio transients in general and core-collapse RSNe in particular 



(IGal-Yam et al.ll2006l ). The unprecedented sensitivity of SKA could allow detection of core- 
collapse RSNe out to a redshift z ~ 5. These detections will be unbiased and automatic in 
that they can occur anywhere in the large SKA field of view without need for targeting based 
on prior detection at other wavelengths. With SKA, the core-collapse RSN inventory should 
increase from the current number of several dozen to ~ 620 yr _1 deg -2 . EVLA should detect 
~ 160 yr -1 deg -2 , and other SKA precursors should reap similarly large RSN harvests. In 
contrast, intrinsically dim RSNe such as Type la events and 1987A-like core-collapse explo- 
sions are unlikely to be found blindly. However, the SKA (and precursor) sensitivities will 
offer the possibility of detecting these events via targeted followup of discoveries by optical 
synoptic surveys such as LSST. 

The science payoff of large-scale RSNe searches touches many areas of astrophysics and 
cosmology. We conclude with examples of possible science applications with the new era of 
RSN survey. However, the true potential of untargeted radio search is very likely beyond 
what we mention. 

Non-prompt RSN emission requires the presence of circumstellar matter, so such sur- 
veys will probe this material and the pre-supernova winds producing it. For core-collapse 
supernovae, pre-supernova winds should depend on the metallic i ties of the progenitor stars 



flLeitherer et all Il992l ; Kudritzki fc Pulsl l2000l ; IVink et all l200ll : Mokiem et all I2007T ). and 



should be weaker in metal-poor environments with lower opacities in progenitor atmospheres. 
This effect should lead to correlations between RSN luminosity and host metallicity, as well 
as an evolution of the RSN luminosity function towards lower values at higher redshifts. For 
Type la supernovae, the mass-loss rate from the progenitors depends on the nature of the 
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binary system, i.e., single or double degenerate (jNomoto et al.lll984J ; llben &: Tutukovlll984l ; 
Webbinklll984l ). Radio detection of Type la supernovae will probe the mass and density 
profile of the surrounding environment and hence be valua ble for studying Type la physics 
feck et al.l E995I : JPanagia et aDbood : Ichomiuk etHlboilh . 



Large-scale synoptic RSNe surveys will complement their optical counterparts. While 
optical surveys such as LSST will provide very large supernova statistics at z < 1, ra- 
dio surveys will be crucial for detections at higher redshifts. The nature and evolution of 
dust obscuration presents a major challenge for optical supernova surveys and supernova 
cosmology. Current s tudies suggest dust obs curation increases rapidly with redshift, but un- 
certainties are large. iMannucci et al.l ( 120071 ) estimate that optical surveys may miss ~ 60% 
of core-collapse supernovae and ~ 35% of Type la supernovae at redshift z ~ 2. RSN ob- 
servations, in contrast, are essentially unaffected by dust. Thus, high-redshift supernovae 
could be detected at radio wavelengths but largely missed in counterpart optical searches. 
Comparing supernova detections in both optical and radio will provide a new and indepen- 
dent way to measure dust dependence on redshift. In particular, SKA will be a powerful 
tool to directly detect supernovae in dust-obscured regions at large redshift, and therefore 
offer what may be the only means to study the total supernovae rate, star-formation, and 
dust behavior in these areas. 

Additionally, radio surveys will reveal rare and exotic events. For example, some 
Type Ibc supernovae are linked to long G RBs (Galama et al. 1998; and see reviews in 
Woosley fc Bloomll2006l ; iGehrels et al.ll2009f ). probably via highly relativistic jets po wered by 
centr a l engines and will manifest themselves in ex t remely luminous radio emission ( Woosley. 



19931 : llwamoto et all Il998l : JLi & Chevalier! Il999l ; IWooslev et al.l Il999l : iHeger et all 12003 ). 



Thus one might expect radio surveys to preferentially detect more Type Ibc supernovae than 
other supernova types. An unbiased sample of Type Ibc RSNe will provide new informa- 
tion about the circumstellar environments of these explosions and thus pr obe the mass-loss 



effects believed to be crucial to the Ibc pre-explosion evolutionary path (|Price et al.l 12002 



Soderberg et alJl2004 12009 : ICrockett et al.ll2007l : IWellons fc Soderberdl201ll ): in addition, a 
large sample of Ibc RSNe will all ow systematic study of the differences, if any, between thos e 
which do an do not host GRBs (jBerger et al.ll2003l ; ISoderberg et al.ll2006t ISoderbergll2007l ). 



Furthermore, radio surveys give unique new insight into a possible class of massive star 
deaths via dire ct collapse into black holes, with powerful neutrino bursts but no electromag 



netic emission (JMacFadyen fc Woosleylll999l ; lFryerlll999l ; iMacFadyen et al.ll200ll ; IHeger et al. 



20031 ). These "invisible collapses" can be probed by comparing supern ovae detected ele ctro- 
magnetically and the diffuse background of cosmic supernova neutrinos (JLien et al.ll20T0L and 
references therein). By revealing dust-enshrouded SNe, radio surveys will make this com- 
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parison robust by removing the degeneracy between truly invisible events and those which 
are simply optically obscured. Indeed, direct collapse events without explosions but with 
relativistic jets are candidates for GRB progenitors. A comparison among RSNe, optical 
supernovae, GRBs, and neutrino observations will provide important clues to the physics of 
visible and invisible collapses, and their relation with GRBs. 

We thus believe that a synoptic survey in radio wavelengths will be crucial in many 
fields of astrophysics, for it will bring the first complete and unbiased RSN sample and 
systematically explore exotic radio transients. SKA will be capable of performing such an 
untargeted survey with its unprecedented sensitivity. Our knowledge of supernovae will thus 
be firmly extended into the radio and to high redshifts. 

We thank Kurt Weiler, Christopher Stockdale, and Shunsaku Horiuchi for encourage- 
ment and illuminating conversations. We are also grateful to Joseph Lazio for helpful com- 
ments that have improved this paper. 
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